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Thermodynamics of gravity favours Weak Censorship Conjecture 
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We use the formulation of thermodynamics of gravity as proposed by Clifton, Ellis and Tavakol 
on the gravitational collapse of dustlike matter, that violates the strong or weak cosmic censorship 
conjecture depending on the initial data. We transparently demonstrate that the gravitational en¬ 
tropy prefers the scenario where the stronger version is violated but the weak censorship conjecture 
is satisfied. This is a novel result, showing the weak cosmic censorship and hence the future asymp¬ 
totically simple structure of spacetime, is being validated by the nature of gravity, without imposing 
any extra constraint on the form of matter. 
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Since the inception of Cosmic Censorship Conjec¬ 
ture (CCC) more than four decades ago [T], there have 
been numerous unsuccessful attempts towards its precise 
mathematical proof, within the context of general rela¬ 
tivity mu. The reason for the continuing endeavour for 
establishing CCC on a solid mathematical footing is ob¬ 
vious; this conjecture is one of the key ingredients for 
the proofs of all important theorems of black hole dy¬ 
namics and thermodynamics [S]. Therefore any trans¬ 
gression of CCC would imply a reformulation of other¬ 
wise well established results on the structure of black hole 
physics, that has led to many remarkable theoretical de¬ 
velopments. Hence comes the necessity of the theoretical 
censoring of naked singularities in nature. 

As of now, there are two versions of CCC that are used 
extensively in the literature [2] . The physical formulation 
of the weaker one (WCC) states that “All singularities of 
a gravitational collapse are hidden within black holes and 
cannot be seen by a distant observer”. The stronger ver¬ 
sion (see) is physically formulated as follows: “All phys¬ 
ically reasonably spacetimes are globally hyperbolic. That 
is, apart from possible initial singularities like Big Bang, 
no singularity is ever visible to any observer”. In the 
context of gravitational collapse of astrophysical compact 
stars, the impact of the above conjectures can be sum¬ 
marised in the following way. The stronger version as¬ 
serts that the spacetime singularity thus formed must be 
immersed in the trapped region and hence there cannot 
be any outgoing future directed non-spacelike geodesics 
from the vicinity of the singularity. Whereas the weaker 
version admits the possibility of the existence of such 
geodesics, but these will definitely enter the trapped re¬ 
gion in future and fall back to the singularity, thus saving 
the future asymptotic structure of the spacetime mani¬ 
fold. 

In spite of the absence of any mathematical proof of 
cosmic censorship, numerous studies have produced a 
large number of counterexamples, where one or both of 


the above versions are violated during dynamical gravi¬ 
tational collapse of different forms of matter fields under 
very generic conditions (see 011 [SI HI for detailed ref¬ 
erences). The main reason for violation of CCC in these 
counterexamples is that the trapped surfaces do not form 
early enough to shield the singularity from outside ob¬ 
servers. One possible way to to avoid the unpleasantries 
of cosmic nudity would be to assert that the form of mat¬ 
ter fields, that lead to the violation of censorship, are 
unphysical. However, these counterexamples span a wide 
class of matter fields (namely perfect fluids, null fluids, 
scalar fields or a combination of these) obeying all phys¬ 
ically reasonable energy conditions and such a strong re¬ 
striction on the thermodynamics of matter helds would 
itself be unphysical. Furthermore, in a recent investiga¬ 
tion 0, it was shown that the Weyl curvature (which is 
a measure of free gravity) plays an important role in the 
violation of censorship. 

Under these circumstances the key question that arises 
here is as follows: Instead of thermodynamics of mat¬ 
ter, should the thermodynamics of the gravitational field 
dictate the final outcome (in terms of a black hole or 
a naked singularity) of the gravitational collapse? Al¬ 
though this is a very important question, until recently 
it would have been impossible to investigate upon, due 
to the unavailability of suitable local definitions of ther¬ 
modynamic quantities of the gravitational field. In fact, 
Penrose 0iini did suggest that the definition should be 
based in properties of the Weyl tensor (instead of the 
Ricci tensor which is directly related to matter helds via 
Einstein held equations) but gave no specihe formula. 
Recently, in a seminal work m, it was proposed that 
the energy density and pressures for the gravitational 
held can be extracted from the square root of the Bel- 
Robinson tensor [T2] , which is constructed from the Weyl 
tensor and it’s dual and acts like an effective energy- 
momentum tensor of the free gravitational held m- It 
was also proposed that the local temperature of the grav- 
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itational field for an observer with 4-velocity can be 
generalised from the expression of the surface gravity |14j , 
which is given as 

Tgrav = bUa ^ (1) 

TT 

where P and are real null vectors in a Newman- 
Penrose tetrad. These definitions have been shown to 
satisfy the necessary conditions in order to generate a 
well defined entropy and have already been applied in 
different scenarios: the analysis of scalar perturbations in 
FLRW spacetime mm, the growth of inhomogeneities 
in perfect fluids, Oppenheimer-Snyder dust collapse and 
its final black hole stage [T^. The results obtained are 
either consistent with previous studies or in line with ex¬ 
pected thermodynamic behaviour m- In particular it 
has been possible to prove that the formation of struc¬ 
tures under gravity is a thermodynamically favoured pro¬ 
cess and that the usual Bekenstein-Hawking entropy is 
the final result of the change in gravitational entropy dur¬ 
ing collapse that eventually forms the black hole. 

Through these definitions, it was shown that it is pos¬ 
sible to define a first law of thermodynamics for the grav¬ 
itational field as 

TgravSgray — il-l-gravV) + PgrayV , (2) 

where Tgrav, Pgrav, Pgrav and Sgrav represent the effec¬ 
tive temperature, energy density, isotropic pressure and 
entropy density of the gravitational field. Here v is the 
spatial volume element and the dot represents the deriva¬ 
tive with respect to time. For Petrov type D spacetimes 
we have Pgrav = 0 and the other thermodynamic quanti¬ 
ties El can be written as 

dS - ^ if^gravV) _ J_|c| /O'! 

^^grav — rp •> h^grav — q 

grav 

where £ is the Weyl curvature scalar constructed via 
the semi-tetrad 1-I-1-I-2 covariant formalism [T5I420) . Fur¬ 
thermore, if we assume spherical symmetry, in the semi¬ 
tetrad formalism the expression of the gravitational tem¬ 
perature Q becomes 


of gravitational entropy between two different configu¬ 
rations SSgrav = Sgravih) - Sgrav{ti), where t is the 
relevant timelike coordinate (an affine parameter to the 
integral curves of u“). 

Now, to clearly see whether the thermodynamics of 
free gravity indeed have a preference towards CCC or 
otherwise, we need an exact model of a collapsing star 
that obeys/disobeys either version of censorship conjec¬ 
ture, depending on the initial data at the time slice from 
which the collapse commences. The simplest example 
is a spherically symmetric collapsing star of marginally 
bound dustlike matter, that can be smoothly matched to 
a vacuum Schwarzschild exterior. In this case, the metric 
in the interior of the star will be the marginally bound 
Lemaitre-Tolman-Bondi metric in comoving (t, r,6,4>) co¬ 
ordinates: 

ds^ = —dt^ + R\t, rYdr"^ -\- R{t, r)^ -I- sin^ dd(j)^) , 

(5) 

where R{t, r) is the area radius of the collapsing spheri¬ 
cal dust shell labeled by r at an instant t and the prime 
denotes derivative with respect to r. By the scaling free¬ 
dom of the radial coordinate, we can always scale the 
coordinate r such that at initial epoch t = ti, from which 
the collapse commences, we have R{ti,r) = r. This col¬ 
lapsing star is smoothly matched to the Schwarzschild 
exterior at the boundary shell denoted by r = rj,. From 
the {t,t) component of Einstein field equations for the 
metric ([^ we get F' = pR^R', where F = F{r) is the 
Misner-Sharp mass [21] of the star which is related to 
the Schwarzschild mass M of the exterior spacetime by 
F(rb) = 2M, and p{t,r) is the density of the collapsing 
star. The mass function can be related to the initial den- 
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sity Pair) = p(ti,r) by the relation F{r) = J po{x)dx. 

0 

Thus we can always write F{r) = r^A4(r), where A4(r) 
is a well defined function Vr G [0, r^j. 

For the marginally bound case, the equation of motion 
for the collapsing shells is given by RR^ = F [551 [53] . 
Solving the equation of motion we obtain 

/ Q _\ 2/3 

R{r,t) =r -y/M(j^tj . (6) 
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where A is the 4-acceleration, 0 is the volume expan¬ 
sion and E is the shear scalar of the congruence of the 
observer’s 4-velocity This quantity can be identi¬ 
fied as the temperature of the gravitational field mea¬ 
sured locally at a point of the spacetime, associated with 
the 2-sphere through that point. Once all quantities are 
specified, one can integrate eq.([^ on a spacelike hyper¬ 
surface to obtain the instantaneous measure of gravita¬ 
tional entropy Sgrav for that specific time slice. In the 
following we will be concerned also with the variation 


To determine the time at which every shell crosses the 
horizon and becomes trapped, from the equation of the 
apparent horizon F = R, we extract the general expres¬ 
sion 


tah{r) 


2 

3y^M{r) 



(7) 


To demonstrate transparently the role of initial data in 
the violation of CCC, let us consider a realistic initial 
density profile: po{r) = poo ~ ■ To avoid patholo¬ 

gies like shell crossing singularities, and to mimic realistic 
stars where the density is a non-increasing function of the 
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radial coordinate, we must have po 2 > 0. This initial data 
will then imply A^(r) = (poo/3) —(/9o2/5)r^ = toq—W 2 r^. 
A central singularity forms at = 2 /while the 
outermost shell rh collapses in the singularity at time 



FIG. 1: The two scenarios depicting globally naked (left) and 
locally naked (right) situations. The red solid line is the out¬ 
going null geodesic and the blue dashed line is the apparent 
horizon in {t, r) plane. 


It has been shown quite rigorously by many authors 
(see for example [11[S1IH] and the references therein) that 
whenever po 2 > 0 the central singularity is locally naked 
with a future directed radial null geodesic emerging from 
it. The equation for this null geodesic tng{r) in the {t,r) 
plane can be obtained by integrating the null geodesic 
equation {dt/dr)ng = R' , with the initial value arbitrarily 
close to the central singularity. The existence of such null 
geodesic implies violation of SCO. However WCC may 
still be satisfied if the null geodesic enters the trapped 
region in the future at r = r*, where r* is the smallest 
real positive root of the equation tng(r) = tahir). Thus, 
with the given initial data, we may have three distinct 
final outcomes of the collapse: 

1. Case I: po 2 = 0. This is the well known case of 
Oppenheimer-Snyder-Datt [211 US] collapse and the 
vicinity of the central shell gets trapped before the 
singularity formation. In this case, there cannot be 
any outgoing null geodesics from the central singu¬ 
larity and hence both SCC and WCC are satisfied. 

2. Case II: po 2 > 0, rt, > r*. In this case the cen¬ 
tral singularity is locally naked, however the future 
directed null geodesic from the singularity enters 
the trapped region sometime in the future and falls 
back to the singularity. In this case SCC is vio¬ 
lated but WCC is still satisfied and the asymptotic 
structure of the spacetime does not change. 

3. Case III: po 2 > 0,rh < r*. In this case the null 
geodesic from the naked central singularity leaves 
the stellar surface before being trapped and trav¬ 
els to future null inhnity, destroying the asymptotic 
structure. In this case both SCC and WCC are vi¬ 
olated and the central singularity is globally naked. 
Figure 1 shows the difference between a locally and 
globally naked singularity. 


To examine which of these three cases are most 
favoured by the thermodynamics of gravity, we calcu¬ 
late the net entropy change from the initial epoch to the 
epoch of the central singularity t = tgg, as measured by 
congruences of free falling observers both in the interior 
of the star and the exterior Schwarzschild spacetime. To 
do this, we must exactly match the time slices of interior 
and exterior spacetimes, which can be done by describing 
the Schwarzschild geometry in Lemaitre coordinates 

= —dt'^ + -I- dfl^ , (8) 

R 

where the area radius is given by R{f, t) = (f (r — 
and the comoving radial coordinate f labels each free 
falling trajectory. 

In order to obtain the gravitational entropy for the ex¬ 
terior spacetime, we need to integrate eq. on a space¬ 
like hypersurface, with gravitational energy density and 
temperature given respectively by 

_ 1 2M _ 1 

where we have removed the absolute values since R and 
M are both positive. Therefore the total gravitational 
entropy over a generic time slice of the exterior spacetime 
reduces to 


'S'IrLW =47rM(foo-rs), (10) 

where fg is the coordinate of the free falling collapsing 
stellar surface, while Poo denotes the free falling observer 
at infinity. This implies that the change in gravitational 
entropy in the exterior spacetime between any two co¬ 
moving time slice as measured by congruences of free 
falling observers vanishes: 


SSt±, = ST±,{t2) - SlfM = 0 - 


grav 


grav v 


( 11 ) 


Now for the interior spacetime we can easily calculate the 
relevant quantities in eq.([^ and eq.Q to get 


1 

l^grav — 7^ 
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( 12 ) 

(13) 


After some manipulations, we can write the gravitational 
entropy at any time slice “t” in the following useful form: 


'b 

5“L(i) = 2^1 


£/p 


f/p + 3 


R' VFr dr. 


(14) 


Figure 2 shows the numerically integrated results of 
<5'S'g”L = Sp^v{iso) - Splviti) as a function of the 
boundary of the star for three pairs of toq and m 2 - 
















4 


We recall that when r* = rb the singularity is marginally 
locally naked. It can be easily seen that the change in 
total entropy is a monotonically increasing function of 
the boundary radius of the collapsing star, for physically 
realistic values of the initial density parameters. 


^Sgrav 



FIG. 2: Numerical integration of SSgravicb) in a neighbour¬ 
hood of r* = 0.2 for different points {mo,m 2 } in the pa¬ 
rameter space: {0.3,0.05} (dotted), {0.2,0.018} (dashed) and 
{0.1,0.003} (solid). 

Now, let us revisit the three distinct cases mentioned 
before, in light of the net gravitational entropy change 
during the process of collapse. We would like to em¬ 
phasise that since dust is an isentropic form of matter, 
there will be no change in matter entropy for the col¬ 
lapsing star. Hence, change in the gravitational entropy 
will account for the total entropic evolution. For the first 
case, where both SCC and WCC are satisfied, since the 
Weyl tensor in the interior spacetime is identically zero, 
there is no change in gravitational entropy, and hence 
this will be the least favoured case. Now to compare be¬ 
tween the second and third cases, we can see from Figure 
2 that the net gravitational entropy change for rt > r* is 
greater than that for rt < r*. Therefore the second case, 
in which SCC is violated but WCC is not, will in general 
be more favoured thermodynamically. 

This is indeed a very interesting result, as it clearly 
shows that the thermodynamics of gravity discards the 
homogeneous case which is clearly unstable under mat¬ 
ter perturbations. Moreover, it is fascinating to see that 
gravity is favouring the scenario that saves the asymp¬ 
totic structure of the spacetime, validating most of the 
important theorems of black hole dynamics and thermo¬ 
dynamics. This result leads us to argue that any general 
mathematical proof of censorship conjectures should nat¬ 
urally emerge from the thermodynamic properties of the 
gravitational held, a detailed study of which will be re¬ 
ported later. 
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